function [a,b,c, absbank] = QuaBankProfile(e, absBank0, absBankmid, absBankf)
a = 2*(absBank0 - 2*absBankmid + absBankf)/(e0 - ef)^2;
b = -(absBank0*e0 + 3*absBank0*ef - 4*absBankmid*e0 - 4*absBankmid*ef + 3*absBankf*e0 + absBankf*ef)/(e0 - ef)^2;
c = (absBank0*e0*ef + absBank0*ef*ef - 4*absBankmid*e0*ef + absBankf*e0*e0 + absBankf*e0*ef)/(e0 - ef)^2;
absbank = a*e^2 + b*e + c;
end